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Abstract 

We analyse g-functional equations arising from tree-like combinatorial structures, 
which are counted by size, internal path length, and certain generalisations thereof. 

The corresponding counting parameters are labelled by a positive integer k. We 
show the existence of a joint limit distribution for these parameters in the limit of 
infinite size, if the size generating function has a square root as dominant singularity. 

The limit distribution coincides with that of integrals of k-th. powers of the standard 
Brownian excursion. Our approach yields a recursion for the moments of the limit 
distribution. It can be used to analyse asymptotic expansions of the moments, and 
it admits an extension to other types of singularity. 

Keywords: simply generated trees, g-difference equation, Brownian excursion, limit 
distribution 


1 Introduction and main results 

1.1 A combinatorial motivation 

When studying combinatorial classes, a functional equation of the form 

G{x) = P{x,G{x)) (1.1) 

frequently arises, where G{x) is the generating function of the class, and P{x, y) is a formal 
power series in two variables with real coefficients. Prominent examples are classes of simply 
generated trees, counted by number of vertices [37] , classes of directed square lattice paths 
counted by length [3|, and classes of square lattice polygons, counted by perimeter [9]. 
Indeed, there exist combinatorial bijections between corresponding models of trees, paths. 
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and polygons, see e.g. [19] for Catalan trees, Dyck paths, and staircase polygons. Eqn. fll.ip 
reflects a combinatorial decomposition of the given class. If P{0,y) = 0 and ^{0,y) ^ 1, 
then eqn. fll.ip admits a nniqne solntion G{x) G lR[[x]] snch that G(0) = 0. Here, M[[x]] 
denotes the ring of formal power series in x with coefficients in M. Often, the coefficients of 
P{x,y) are non-negative real nnmbers. Then, the series G{x) is a power series with non¬ 
negative coefficients, typically analytic at x = 0 with a sqnare root as dominant singnlarity, 
see e.g. m Thm. 10.6], [16], or pTl Ch. VII.4], and references therein. 

We are interested in certain deformations of the above eqnation. This is done by 
introdncing a new variable q, snch that the limit g —> 1 rednces to the original eqnation. 
For example, the fnnctional eqnation 

G{x, q) = P{x, G{qx, q)) (1.2) 

defines a formal power series in x with polynomial coefficients, i.e., G{x,q) G M[g][[x]], 
if eqn. fll.ip defines a power series G{x) G M[[x]]. The above eqnation is a q-dijference 
equation, see e.g. [51] and references therein. It appears in classes of simply generated 
trees, connted by nnmber of vertices and internal path length (i.e., the snm of the vertex 
distances to the root), in classes of directed sqnare lattice paths, connted by length and 
area nnder the path, and in classes of sqnare lattice polygons, connted by perimeter and 
area. For some models, an explicit expression for its generating fnnction G{x, q) is known, 
see e.g. [91111 US] and references therein. Snch an expression typically contains g-prodncts 
and has a natnral bonndary |g| = 1. An interesting qnestion concerns the statistics of the 
additional connting parameter, e.g., in a nniform ensemble in the limit of large system size. 
It is known ([50l |19], see also [16]) that eqn. (11.2p leads, for certain g-difference eqnations 
and after appropriate normalisation, to the Airy distribution as the limit distribntion for 
the additional parameter. This distribntion is known to also describe the area nnder a 
Brownian excnrsion, see the following snbsection. Note that generalisations of eqn. fll.ip 
other than eqn. fll.2l) have also been stndied previonsly. A class of eqnations, which leads 
to Ganssian limit laws, is discnssed in [TB] . 

The idea of iterating the above deformations has been considered by Dnchon [T8l [T9] . 
The deformation variables may be denoted by g^, where fc G {1,..., M}, and the associated 
connting parameters are called parameters of rank k + 1. For eqn. fll.2p . an example is 
given by 

G{x, gi,..., qn) = P{x, G{xqi ■ ...-qM, gig 2 • ... ■ gvf, q 2 q 3 • ... ■ gM, • • •, qn))- (1-3) 

Here, G{x, gi,..., qu) = '^nPni.Qiy ■ ■ ■ y Qm)x'^ is a formal power series in x, with polynomial 
coefficients Pn{qi, ■ ■ ■, qn) £ ^[qi, ■ ■ ■, Qm]- The name “parameter of rank k + 1” reflects 
that k is the smallest integer r, snch that the degree of the polynomial Pn(?i, ■ ■ ■, qu) in 
qk is bonnded by for some constant c (see [T9l Lemma 1] and [H]). We will call 

an eqnation like eqn. fll.Sp a q-functional equation fDefinition 12.41) . Again the qnestion 
arises, nnder which conditions a limit distribntion for the additional connting parameters 
exists. In this paper, we shall show this for a class of deformations which we call g- 
shift fDefinition 12.11) . the main assnmption on the g-fnnctional eqnation being that the 
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solution of the undeformed equation, eqn. fll.ll) . is analytic at the origin, with a square 
root as dominant singularity fAssumption 14.ID . See Theorem 11.51 for a precise statement. 
The resulting limit distributions appear to be related to distributions of integrals of fc-th 
powers of the standard Brownian excursion. We will obtain a recursion for the moments 
of the joint distribution. Our approach is based on the multivariate moment method (see 
e.g. [271 d). The univariate case M = 1 has been studied previously [501 IMl HSl SS] , and 
recursions for M = 2 have been studied in [iQ]. 

Before we consider general g-functional equations in Section [2l let us hrst discuss the 
above questions in more detail for the simple example of Dyck paths. 

1.2 Dyck paths and Brownian excursions 

We review the connection between Dyck paths and Brownian excursions. This relates, 
in particular, the moments of height of a random Dyck path to Brownian excursion mo¬ 
ments. Since the generating function of Dyck paths, counted by length and moments of 
their height, provides a simple example of a solution of a g-functional equation, this also 
motivates the appearance of excursion moments in g-functional equations. We will state 
a central result of the paper in Theorem 11.21 which characterises the Brownian excursion 
random variables by a recursion for their moments. 

Let No = N U {0} and M>o = {a: G M : a: > 0}. A Dyck path of length 2n, where n G Nq, 
is a map y : [0,2n] —>• M>o, where j/(0) = y{2n) = 0 and \y{i) — ?/(« -|- 1)| = 1 for z G No 
and i < 2n. For non-integral argument, y{s) is dehned by linear interpolation. The values 
y{s) are called the heights of the path. An arch of length 2n is a Dyck path of length 2n, 
where n > 0 and y{s) > 0 for s G (0, 2n). An example is given in Figure [H 



Figure 1: A Dyck path of length 2n = 10. It is a sequence of two arches of lengths 8 and 2. 

We are interested in a probabilistic interpretation of Dyck paths in a uniform ensemble 
where, for hxed length, each of the hnitely many Dyck paths occurs with the same proba¬ 
bility. For 0 < s < 2n, let W(s) denote the height of a random Dyck path of length 2n. It 
is well known, see e.g. [13], that the average maximal height h{n) of a random Dyck path 
of length 2n has the asymptotic form h{n) ~ -\/Fn as n —> cxo. In order to obtain a hnite 
positive limit as n approaches inhnity, and in order to normalise the domain, we introduce 
the normalised height Ynit) = (2n)“^/^W(2nt), where 0 < f < 1. 
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The sequence {F„(t)}neN is a sequence of stochastic processes defined on ((^[0,1], 11 • | |oo) 
with the Borel cr-algebra, which converges in distribution to the standard Brownian excur¬ 
sion e{t) of duration 1, see [6],[3l[36]. This implies convergence in distribution of sequences 
of continuous bounded functionals of Ynit) towards the corresponding excursion function¬ 
als. (We refer to [8] for background about convergence of probability measures.) There is 
a more recent result m Thm. 9] that also 

hm E„[F(W(t))] =E[F(e(f))] 

n—>oo 

for continuous functionals F : C[0,1] —>■ M of polynomial growth, i.e., for functionals 
such that there exists an r > 0 with \F{y)\ < ||?/||^ for all y G ^[0,1]. The above 
property is called polynomial convergence. In particular, polynomial convergence implies 
moment convergence for such a functional, i.e., convergence of the sequence of moments 
of order /, {E„[F(W(f))*]}„gi^, for every I G Nq. For the functional of polynomial growth 
F{y) = y^{t)dt, fc G M, we call the random variable e’^{t)dt the k-th excursion 
moment. 

As counting parameters for Dyck paths of length 2n, we consider the parameters 

2n 

Xk,n = '^y''{i) {k = l,...,M). (1.4) 

i=0 

These are sums of fc-th powers of heights, which we call k-th moments of height. The 
parameter is the area under the Dyck path, the parameter X 2 ^n is called the moment 
of inertia of the Dyck path [10] • Let denote the fc-th moment of height of a random 
Dyck path of length 2n. In terms of W(s), the random variables Xk^n are expressed as 

2n 

i=0 

The scaling of the average height of a random Dyck path with its length suggests consid¬ 
ering normalised random variables dehned by 


( Y Y Y i — I ^2,n Xm^ti \ 

^ 2 , n , ■ ■ • , -^M,n) — I ( i + 2)/2 ’ „( 2 + 2)/2 ’ ' ' ' ’ „( M + 2)/2 ) ' 


(1.5) 


Then, in terms of Yn{t), the normalised random variables Xk^n are given by 

2n 

Y’^l, 

2n j 2n 


Xkn = 2(^+2)/2 ( J_y 


i=0 


( 1 . 6 ) 


Since the above sum is a Riemann sum, one is led to expect convergence in distribution 
and moment convergence of the sequence to e^{t)dt, due to the con¬ 

vergence properties of {W,(t)}„gis}. This is the statement of the following proposition. 
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Proposition 1.1. For fc G {1,..., M}, let Xk = fg e^(t) dt denote the k-th excursion mo¬ 
ment. The sequence of Dyck path random variables eqn. fll.Sp converges to the normalised 
excursion moments ..., in distribution, 

X 2 ,„, . . . , XM,n) ^ ( 2 ( 1 + 2 )/ 2 ^^^ 2(2+2)/2X2, . .., 2^^+^^/^Xm) (n ^ CX)). (1.7) 


We also have moment convergence. 

Proof. For k G {1,...,M}, define Yk^n = fo^n(i)dt. As stated above, convergence 
in distribntion and moment convergence holds for the seqnence {(Yi 
see IS i and m Thm. 9]. We argne that the seqnences {(XiXM,n)}neN and 
|( 2 (i+ 2 )/ 2 y-^^^^... ^ Converge in distribntion and for moments to the same 

limit. 

Fix k G {1,...,M} and consider, for a Dyck path y of length 2n and for m G 
{0,1,..., 2n — 1}, the elementary estimate (max{0, y{m) — 1})^ < y^{s)ds < {yijn) + 
1)^. Snmming over m yields, together with the binomial theorem, the estimate 

E £ f s E (1-8) 

This estimate translates to the distribntion fnnctions of the corresponding random vari¬ 
ables. In terms of the random variables Xk^n and Yn(t), we get 


Z=0 ^ 2 -'U z=o ^ 2 


X, 


Ln 


nik-l)/2 ’ 


(1.9) 


where the snperscript indicates that the relation is to be nnderstood via distribntion fnnc¬ 
tions. Since the seqnence of random variables {Xi^n/n!'^~^'’^‘^}nm converges in probability 
to zero, for I = 0,..., k — 1, the random variables on the l.h.s. and on the r.h.s. of eqn. fll.9p 
converge in distribntion to the same limit. We conclnde that the seqnences {Xfc^„}„gN and 
| 2 (fc+ 2)/2 Y^(t)dt}„gN converge in distribntion to the same limit. Moment convergence 
follows similarly from eqn. fll.Sp . 

Since the above estimates also hold jointly in /c G {1,... ,M}, the statement of the 
proposition follows. □ 

Remark. The above statement concerns Dyck paths in a nniform ensemble. The same 
argnment as above leads to analogous results for more general ensembles of Dyck paths, 
which arise from the depth hrst process of simply generated trees Compare Theo¬ 

rem 11.51 for a further generalisation. 

Due to the above result. Brownian excursion functionals may be studied via their 
discrete Dyck path counterparts. Below, we will analyse the asymptotic behaviour of 
the moments of the random variables (Xi „,... ,XM,n) eqn. fll.5|] . using the ^-functional 
equation which the generating function of Dyck paths obeys. This implies a certain re¬ 
cursion for the moments of the joint distribution of the excursion moments (Xi,.. . , Xm)- 
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We have the following result, which is obtained by combining the statements of Proposi¬ 
tion 11.11 and Proposition 11.41 below. For M-dimensional vectors, we use the abbreviation 
k = {ki ,..., /cm) and write 0 = (0,..., 0). For i G {1,..., M}, the unit vector e* is dehned 
by {ei)j = 6ij for j = 1,..., M, and Z < fc for vectors Z = (Zi,..., Im) and k = (fci,..., Zcm) 
means that h < ki for i = 1,..., M. A multivariate power series with complex coefficients 
is called entire, if it converges for arbitrary complex arguments. 


Theorem 1.2 (Excursion moments). For fc G {1,..., M}, let e^(t) dZ denote the 

k-th excursion moment. The moments of the joint distribution of (Xi,... ,Xm) are given 
by 


■ ... ■ W] _ ^ . . 

ki\-...-kM\ r(7fci,...,fcM) 2 


( 1 . 10 ) 


where 'yki,...,kM = 7fc = ~l/2 + + V2)^i; and where F( 2 ;) is the Gamma function. 

For k 0, the numbers fki,...,kM = fk are characterised by the recursion 


7fc—ei fk 


-ei 


M-1 


2(* l){ki + l)fk- E fifk-i — 0, 

i=l 0<l<k 


( 1 . 11 ) 


with boundary conditions fo = —4 and fk = 0, if kj < 0 for some j G {!,..., M}. The 
moment generating function is entire. Hence, the joint distribution is 

uniguely defined by its moments. □ 

Remarks, i) The theorem asserts that the numbers fk can be recursively computed from 
eqn. fll.lll) . To see this, we remark that the above equation has — fk = 2fkfo/S as a term, 
and that all other numbers fi in the equation satisfy I -< k for a suitable total order 
which is specihed below in Dehnition 12.21 

a) In probability theory, Louchard’s theorem [3l] leads to a characterisation of a certain 
Laplace transform of the moment generating function, for some excursion functionals in¬ 
cluding excursion moments. It is however difficult to extract moment values or moment 
recursions from these expressions. For M = 1, this has been done in [55]. The moment 
generating function of the marginal distribution for M = 2 has been obtained in [401 
Thm. 2.4]. Moment values or recursions for M > 2 have apparently not been derived. 

Hi) We obtain Theorem 11.21 bv studying the corresponding Dyck path functionals. Our 
discrete approach has been used previously [10|. It led to a combinatorial derivation [40], 
Thm. 3.1] of Louchard’s formula for integrals over Brownian excursion polynomials. It 
also led to the values E[(Xi)^i(X2)^2], see [40] Table 2]. For arbitrary M, our result is 
announced in [47] . [A factor 1/2 is missing on the r.h.s. of eqn. (6) in |47j]. 
iv) Assuming only convergence in distribution of (Xi^„,..., XM,n), the above result can be 
used to provide an alternative proof of moment convergence. 

Let us briefly discuss the moments of the marginal distributions for M < 4. For 
M = 1, the hrst few coefficients E[(Xi)*’] are 1, These 

are related to the Airy distribution [aiianiiaii, i.e., the distribution of area under the 


6 















Brownian excursion \/8e{t). Explicit expressions are known for the moment generating 
function, for the density of the distribution, and for the moments. 

For M = 2, the Erst few moments E[(X 2 )^] are 1, i, i||§. An 

explicit expression for the corresponding moment generating function 001 is 


E|e“=l = X; 


E|(X2 


k=0 


k\ 


k^ = 


(v^)) 




3/2 


The hrst few moments E[(X 3 )*^] are 1, 
the hrst few moments E[(X 4 )^] we have 1, 


3V^ 

16 


207 11907v^ 88655283 1165359069^271 


2§i 


?fll® ’l9?)W°°10516W2'0^3°°® 


2’ 840’ 1201200’ 76236160’ 3259095840000 

an open problem to End explicit expressions for the moment generating functions 


It remains 


1.3 Dyck paths and g-functional equations 

We discuss the functional equation which the generating function of Dyck paths obeys, 
when counted by length and fc-th moments of height, k = 1,..., M. We then state in 
Theorem 11.51 a convergence result for the limit distribution of counting parameters related 
to general g-functional equations. This is the main result of our paper. 

Let u = (uq, Ml,..., um) denote formal (commutative) variables. For a Dyck path p of 


length 2 m, its weight w{p) is given by w{p) = ■ ... ■ where Xk,n is the fc-th 

moment of height eqn. fll.4|) . For the set T) of Dyck paths, its generating function is the 
formal power series D{u) = 'Yhp^v'^iP)- arches, its generating function 

is the formal power series A{u) = variables Mi, ... ,mm are interpreted as 

deformation variables. They may all be set equal to unity, in which case the generating 
function reduces to that of Dyck paths (arches), counted by length only. The generating 
functions satisfy the following functional equation. 

Theorem 1.3 (Dyck path generating function, cf. [ID])- The generating functions of Dyck 
paths D{u) and of arches A{u) satisfy 


D{uq, ..., Um) — 1 + D{uo, ..., um)A{uo, ..., mm), 
A{Uo, Ml, ... , Um) = UqUi ■ ... • umD{vo, Vm), 


( 1 . 12 ) 


where the monomials Vk = Vk{u) are given by 

M 

Mo(n) = M0M1M2 ■ . . . ■ M^, t;;.('u)=JJ 


U 


(fc = l,...,M). 


l=k 


□ 


Remarks, i) For M = 2, the above theorem appears in [IHl Sec. 3.1]. For arbitrary M, 
it is (incorrectly) stated in [IHl P- 713]. The proof for M = 2 given in [ID] generalises 
to arbitrary M. It uses a last passage decomposition argument and the additivity of the 
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counting parameters with respect to the sequence construction, yielding the hrst equation 
in the theorem. Note that the arch decomposition of Dyck paths, as depicted in Figure 
[H implies the equivalent equation D{u) = 1/(1 — A{u)). For n G N, there is a bijection 
between the set of arches of length 2n and the set of Dyck paths of length 2n — 2, by 
identifying an arch, where the bottom layer has been removed, with the corresponding 
Dyck path. This implies, after a short calculation (see also ^81 Thm. 1]), the second 
equation in Theorem 11.31 

a) The induced functional equation for the generating function E{u) = D{u) — 1 is an 
example of a g-functional equation as in Dehnition l2.4l with a square-root singularity in the 
generating function in the “undeformed” case (ui,..., um) = (1, • • •, !)• See the following 
section and Section [HI 

For a probabilistic interpretation of the counting parameters, as in the previous subsec¬ 
tion, let Pno,ni,...,nM denote the number of Dyck paths of length 2nQ, where the fc-th moment 
of height has the value n^, for fc = 1,..., M. We clearly have 0 < 

for no > 0. In the uniform ensemble, the random variables which assign the fc-th 

moment of height to a random Dyck path of length 2no, have the joint distribution 


lP(^l,no ^1) • • • ) ^M,nQ ^m) 


Pno 


Pno,ni,...,nM 

The moments h^fc^^...^fcJ^^(no) of the joint distribution are given by 


(1.13) 




V 


riM 


ki 

rii 


■ IT'M Pno,ni,...,nM 


In the previous subsection, we introduced the normalised random variable eqn. (II.bh . 


/ Y Y Y i — I ^^’^0 ^2, no ^M,no 

Wl,noi ^2,no5 ■ • • 1 — I (i+2)/2’ (2+2)/2 ’ ‘ ‘ ‘ ’ (M+2)/2 

\nQ Hq Hq 

The normalised random variable (Wi^„g, ... ,XM,no) has moments given by 

/ \ '^fcl,...,fcM (^o) 

— (l+2)fci/2+(2+2)A:2/2+... + (M+2)fcM/2' 

Uq 

We argued above that these numbers should tend to a hnite limit, as uq approaches inhnity. 
As we will see below, this is indeed the case. Hence we may dehne 



mku...,kM = hm mfci,...,fcji^(no). (1.15) 

no—^■oo 

A careful analysis of the functional equation eqn. fll.121) . which will be performed in the 
general case from Section jH onwards, yields the following result for the numbers mk-^^,,,^kM- 
Its proof is deferred until Section [HI 
















Proposition 1.4. The normalised moments of Dyck paths eqn. fll.lSp are given 

by 

'OOki,...,kM _ 1 r(^o) , . 

where Mc = 1/4, where r(z) denotes the Gamma function, and where the numbers 7 ^ and 
fk are defined in eqn. (11.111) . The moments have an entire exponential generating function. 
Hence, they define a unique random variable with moments mk. 

Remarks, i) Proposition 11.41 implies convergence in distribution and moment conver¬ 
gence of the normalised random variables eqn. (11.141) . This yields an alternative proof of 
the convergence statement of Proposition 11.11 However, it does not establish a connection 
between the limit random variable and the excursion moments, as in Proposition 11.11 
a) Since the moments mk dehne a unique random variable, the result for the excursion 
moments eqn. fll.lOp follows from Proposition 11.41 bv Proposition ll.il 


For the solution of a general g-functional equation (Dehnition l2.4p . we obtain a similar 
result. Write G{u) = Pno,ni,...,nM'^o°'^T ■ • • • ■ for such a solution. Under mild 

assumptions, the corresponding random variables in eqn. fll.lSp are well-dehned, and the 
same asymptotic analysis as above can be performed. We have the following theorem which 
is, in conjunction with Theorem 11.21 the main result of our paper. Its proof is deferred 
until Section [7] and Section [HI 

Theorem 1.5 (Limit distribution for g-functional equations). Let a q-functional equation 
Definition with solution 

G{Uq, Ml, . . . , Mat) ^ ^ Pno,ni,...,nM^O ^1 ' ■ ■ ■ ' Uj^ 

as in Proposition \2.5\ be given, and let Assumption \4.1\ be satisfied. Assume that the numbers 
Ai in Proposition \4 . 5\ are positive, i.e., > 0 for i = 0,... ,M — 1. Then, the random 

variables (Xi^„q, ..., XM,no) oqn. fll.lSp are well-defined for almost all Hq. The following 
conclusions hold. 

i) The numbers 111-151) . which are derived from the moments of the nor¬ 
malised random variables ..., XM,no) fll.l4p . define a unique random vari¬ 

able (Yi,..., Ym) with moments mki,...,kM- have convergence in distribution, 

(7^1,no5 • • • ! XM,no) (Ih, ■ ■ ■, Ym) {no —>■ cxd), 

and we have moment convergence, 
a) The limiting random variable (Yi,..., Ym) is explicitly given by 

(Yi,..., Yfc) = (ciXi,..., cmXm), 
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where the constants Ck > 0 and the random variables Xk are 



I 


1 


e\t)dt = 


Here, the numbers fVi > 0 are defined in Proposition \4-5[ and e(t) denotes a standard 
Brownian excursion of duration 1. 

Remarks, i) A recursion for the moments of the limit distribution appears in Theorem II .21 
above, see also Proposition 14.51 below. 

ii) For Dyck paths counted by length and fc-th moments of height, the above statements 
follow already from Proposition ll.il More generally, as was argued in the remark following 
Proposition ll.il the above result can be shown to hold for models of simply generated trees 
[37] , counted by number of vertices and fc-th moments of internal path length. This follows 
from the polynomial convergence of the depth hrst process derived from simply generated 
trees towards the Brownian excursion [SI El ESI E|. For a general ^-functional equation, 
such a connection is not known to exist. 

Hi) Our method of proof also allows to study corrections to the asymptotic behaviour, 
compare the discussion in Section [5| and [IS] for the case M = 1. These cannot be 
obtained by the methods of [36] . 

1.4 Structure of the paper 

The remainder of the paper is organised as follows. In Section [21 we introduce g-shifts 
iDehnition 12. ip and g-functional equations iDehnition 12.4p . Our results rely on an appli¬ 
cation of the multivariate moment method (see e.g. [271 Ch. 6.1.] or [71 Sec. 30]). Hence 
in Section El eqn. fl2.9p . we introduce factorial moment generating functions, which are 
derivatives of the solution of the functional equation, evaluated at ui = ... = um = 1- We 
study their properties by a combinatorial analysis of derivatives of the functional equation, 
using a multivariate generalisation of Faa di Bruno’s formula |T5|. In Section |31 we study 
the singular behaviour of the factorial moment generating functions, in the case of a square 
root singularity as the dominant singularity of the size generating function. Proposition 14.51 
gives a recursion for the amplitudes, which describe the leading singular behaviour of the 
factorial moment generating functions. Our method is also called moment pumping [23] . 
We employ in Section [5] an alternative (rigorous) method to obtain the recursion, which 
originates from the method of dominant balance of statistical mechanics [HI EBj. This 
method is generally easier to apply than an analysis of the functional equation via Faa di 
Bruno’s formula, and yields an algorithm for obtaining the amplitude recursion. In addi¬ 
tion, the method allows to analyse corrections to the limiting behaviour. The behaviour 
of the moments follows then by standard methods from singularity analysis of generating 
functions [22l [2l] . In Proposition 15.41 we give a quasi-linear partial differential equation 
for the generating function of the amplitudes. Growth estimates for the amplitudes (and 
hence for the moments) are obtained in Section El by an analysis of the (singular) partial 
differential equation. Existence and uniqueness of a limit distribution is then guaranted by 
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Levy’s continuity theorem, see Section [71 The connection to Brownian excursions follows 
in Section [HI by a comparison with Dyck paths. Possible applications of our method are 
discussed in a concluding section. 

2 Q'-functional equations 

Let C[[tt]] denote the ring of formal power series with complex coefficients in the (commuta¬ 
tive) variables u = (mq, Mi, ..., um)- Let Clw] denote the ring of polynomials with complex 
coefficients in the variables u. We set Uq = (mq, 1, • • •, !)• For u = (mq, Mi, ..., um) and 
n G we dehne to be the monomial = Uq° ■■ u^. We employ the notation 

u+ = (ui,... ,um), the plus sign indicating that the hrst component of u is omitted. For 
n G N and r > 0 , let denote the open polydisc 

= {(xi,..., Xn) G C" : \xk\ < r for all /c = 1,..., n}. 

Let Hr{x) denote the ring of power series in a; = (xi,... ,x„) with complex coefficients, 
which are convergent in D”. 

Definition 2.1 (g-shift). Fix M G N. Let for k = 0,..., M formal power series Vk{u) G 
Clin]] be given, and write v = v{u) = {vq{u),vi{u), ... ,vm{xl)). Assume that there is a 
number d = d{y) satisfying 1 < d < oo, such that 

vq{u) eHd{ui,... ,um)[[uo\], Vk{u) eHd{uk,... ,um) (fc = l,...,M). 

Assume that v{u) satisfies 

dv 

v{uo)=Uq, vo{0,u+) = 0, ^(no) = 1 (fc = 0 ,...,M). 

OUk 

Let r = r{v) be a number 0 < r < cx) such that 

{(ni(n+),...,nM(?i+)) : u+ e Df} C . 

Then, v is called a g-shift, d{v) is called the domain of v, and r{v) is called the range of 

V. 

Remarks, i) For M = 1, a simple example of g-shift appears in the g-difference equation 
eqn. (11.21) . with vo{uo,ui) = uqUi and ni(-ui) = ui. This motivates the name for the 
generalisation. 

a) For a g-shift v, its k-th component Vk{u) does not depend on ui, where I = Q,... ,k — 1, 
and it does depend on Uk. Since v{uo) = Uq, one may interpret v(u) for n 7 ^ no as a 
“deformation” of Uq. The condition no(0,n+) = 0 is imposed to ensure that composition 
of formal power series is well-dehned, see below. 

iiij The identity id : n 1 —^ n is a g-shift with inhnite domain and range. For two g-shifts v 
and w, their composition v ow is well-dehned if r(w) < d(v). As is readily checked, vow 
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is a g-shift in that case, with domain d{v o w) = d{w) and range r{v o w) = r{v). If for a 
g-shift V we have r{v) < d{v), we can thus consider iterated g-shifts where 

= id, = V o (n G M). 

The power series is a g-shift for n G Nq. 

iv) An important subclass (with infinite domain and range) are polynomial g-shifts, i.e., 
g-shifts V satisfying Vk{u) G C)^] for k = 0,, M. Examples are given by the monomial 
g-shifts 

Vk{u) = 

where G = 1 and (n^); = 0 for Z = 0,..., fc — 1 and fc = 0,..., M. These 

appear in eqns. fll.2p and fll.Sp . and for Dyck paths in eqn. fll.121) . In these examples, we 
have (n;);_|_i 7 ^ 0 for Z = 0 ,..., M — 1 . 

For a formal power series G{u) G 7-Zd('u+) [[mq]] ^ *C[[it]] and a g-shift v satisfying 
r{v) < d, we define H{u) G 7-Zrf(^,)(^^+) [[mq]] by 

H{u) =G{v{u)). 

This is well-defined since Hd{v^{u)) C T-Cd{v){u^) and C[[no('ii)]] ^ TZ(i(^,)(^^+)[[^^o]], due 
to ^ 0 ( 0 , 14 +) = 0. We are interested in derivatives of H{u). For clarity of presentation, 
we will use the multi-index notation, and Greek indices /x, u, p will denote vectors with 
non-negative integer entries. For F(u) G C[[ix]] and ly = {uq, ..., pm) ^ we write 

the derivative Fj^{u) G C[[rx]] of F{u) of order ly as 

where we use the convention that Fq{u) = F{u). If G{u) G 7d(i(n+)[[Mo]] C C))^]], we also 
have Gi,{u) G 7-Zd('u+)[[Mo]]- Set \iy\ = pq + ... + pm and iy\ = pq\ ■ ... ■ pm^-- We dehne a 
total order in as follows. 

Definition 2.2 (total order -<). For p = (/xq, ..., Pm) £ Ng"'"^ and v = (z/q, ..., pm) G 
Mg"'"^, we write p ^ u if either |^i| < \iy\, or if |/x| = \iy\, there exists an index k G 
{0,..., M}, such that pk > Pk and pi = Pi for i = 0,..., k — 1. 

The total order introduced above will be used to label terms appearing in derivatives 
of H{u). We have the following lemma. For i G {0,..., M}, let e* G denote the unit 

vector in direction i, given by (ej)^ = Si^k for Zc = 0,..., M. 

Lemma 2.3. Let G{u) G 7-Zd(n+)[[Mo]]; and let a q-shift v with domain d{v) and range 
r{v) < d be given. Set H{u) = G{v{u)) G 7-Zd(i,)('it+)[[Mo]] and fix ly 0. 

i) We have H^{u) G Hd(v){'tJ'+)[[uo]]. For every p satisfying 0 ^ p ^ ly, there exists a 
coefficient A{p,u) G Hd(v){u+)[[uo]], independent of the choice ofG{u), such that 

Hu{.u) = ^ G^j,{v{u)) ■ A{p,u). (2.1) 
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ii) The coefficient A{fj,,u) in eqn. fl2.1l) is, for fi = u, given by 


■4('^.“) = n(sS;(“)) ■ (2-2) 

In particular, we have A{u, Uq) = 1. 

Hi) Fix i G {0,..., M — 1}. If i/j+i > 0, we have u — Cj+i + e* -< i/. The coefficient 
A{pi, u) in eqn. fl2.1l) is, for fx = u — ej+i + e, and u = Uq, given by 

. dvi , . 

A[u - Bi+I + Cj, Uq) = h>i+i- -(no). 

a-Uj+i 

iv) For real numbers r, ro,..., rM, where r^+i > r^ for k = 0,..., M — 1, define = 
^ + J2k=o'^kHk- Then, for indices pL u in eqn. (12.Ih satisfying A{pL,u) ^ 0, we 
have Ofj, < a^. 

Remarks, i) The coefficient A{pi, u) in eqn. fl2.ip might be chosen to vanish. E.g., £x n = 
(1,1) and consider pi = (0,1). We have pi ^ ix but might choose A{pi, u) = 0, as is readily 
verihed by an explicit calculation, using the g-shift property Vk{u) G C[[uk, ■ ■ ■ ,um]]- 
ii) The property iv) will be used for exponent estimates in Proposition 14.41 below. 

Proof, i) This is an application of the chain rule. Successive differentiation leads to 

Hu{u)= ^ G^,{v{u)) ■ A{^i,u), (2.3) 

o<lnl<|i^| 

where A{pL, u) is independent of G{u). To analyse the effect of the g-shift property Vk{u) G 
'kid(v){uk ,..., Um)) consider for fc G {0,..., M} the equation 



It shows that derivatives of II{u) w.r.t. do not contribute to derivatives of G{y{u)) 
w.r.t. Ufc+i,..., vm- a similar statement holds for higher derivatives. This implies that the 
numbers Vk in eqn. fl2.3p can contribute to the numbers fiQ,..., p,k only. If \pi\ = |i/|, we 
thus have < /tq + • • • + for fc G {0,..., M}. We show that pi F v. Assume w.l.o.g. 
that |^i| = li^l. If po > we have pi -< u, and the claim follows. Otherwise, we have 
/To = no. Thus, Vk does not contribute to /xq for fc G {!,..., M}. The previous argument 
yields that z/^ < pi + ... + for fc G {!,..., M}. Now repeat the above argument until 
fiQ = ... = fiM-i = k'o = ... = vm- 1 - Then since |^| = \u\. Thus pi = u, and 

the statement is shown. 

ii) This is seen by an explicit calculation using eqn. fl2.4p . Successively applying derivatives 
w.r.t. Mo; Wi, ..., um yields eqn. fl2.2p . Together with the g-shift property dvk/dukiuo) = 1, 
it follows that A{h',UQ) = 1. Statement Hi) is shown by an analogous calculation. 
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iv) Consider first the case |^| = \u\. Note that, by differentiation, the nnmbers Vk do 
not contribnte to /i^+i,..., /tm, ioi k E {0,..., M}. Thus, the contribution of Uk to is 
maximal if fik = for = 0,..., M. We get ^ = u and = a,^. Since this maximum 
is unique, n -< u implies < a^. Let now |/x| < \u\. For /c G {0,..., M}, denote by Z 4 
the number of derivatives w.r.t. ma,, which contribute to /x. Set v = (z/q, ... ,Fm). Then 
clearly |/x| = \u\ and ap < a^,. The reasoning for the case |^t| = \i'\ can now be applied 
to the present case, with u replaced by v. This yields Thus and the 

statement is shown. □ 

Definition 2.4 (^-functional equation). Let P{u,yi,... , 1 /^) be a formal power series 
P{u, y) G 7-frf(n+) [[uq, y]\, for a number d satisfying 1 < d < oo. Assume that P{0, 0) = 

0, and that §^{0, 0) = 0 for j = 1,..., iV. Let a q-shift with domain d{vA)) 

and range r{v^A'j < j = 1,..., iV. 

i) If the above assumptions are satisfied, the equation 

G{u) = P{u, W'\u ),.... //<">(>*)). ( 2 . 8 ) 

where = G{v^A(^u'^'^ Jqj- j = , iV, is called a g-functional equation. 

a) Let the above assumptions be satisfied. If there is a number d{G) such that 

1 < d{G) < min 

and a formal power series G{u) G 'Hd{G){u+)[[uo]] satisfying eqn. (12.51) . then G{u) is 
called a solution of the q-functional equation. 


Remarks, i) An example of a g-functional equation for M = 1 is given by the g-difference 
equation eqn. fll.2l) . This motivates the name for the generalisation. Examples for M > 1 
appear in eqn. fll.3l) . Examples of g-functional equations frequently satisfy P{u,y) G 
C[n, y] with d = oo, and with monomial g-shifts. We infer from the functional equation 
eqn. fl 1.121) for Dyck paths that the power series E{u) = D{u) — 1 satishes the g-functional 
equation 

E{u) = uoui ■... ■ um iE{u) + 1) {E{v{u)) + 1). (2.6) 

Specialising to u = Uq yields a quadratic equation, which can be explicitly solved for 
E{uq). We get the well-known result 


E{uo) 


1 - 2uo - a/I - 4mo 
2uo 


a) If Q{u,y) G C[[m,j/]] satishes Q(0,0) = 0 and ^{f),y) ^ 1, a formal power series 
G{u) G C[[n]] such that G'(O) = 0 is uniquely dehned as the solution of the equation 
G{u) = Q{u,G{u)). When studying such equations, we may assume without loss of gen¬ 
erality that ^(0,?/) = 0. The above dehnition reduces to this setup, when restricting to 
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U = -Uq. 


A result about solutions of g-functional equations is given by the following proposition. 
We use the vector notation H{u) = , H^^\u)). 


Proposition 2.5. The q-functional equation of Definition 2.4 


G(u)=P(u,H(n)), 


(2.7) 


has a unique solution G{u) G TCd(p){u^)[[uo\] satisfying G(0,n+) = 0, where d{P) = 

Proof If G{u) = T.no>iPnoi'^+)'^o°, then = Eno>i°- 

For no > 1 fixed, we take the coefficient of no“ in eqn. fl2.7p . Due to nh'>(^0,w+) = 0 for 
j = 1 ,..., iV and the assumptions on the derivatives of P{y, u), we get the expression 

Pno{u+) = Wno (^U+, |pi(n5V))> • • • , ( 2 - 8 ) 

for a power series Wn^ {u+,p) G TLd(p){ujf) [p] in M + N{nQ — 1) variables. Thus Pnoiujf) 
is determined recursively in terms of pi{ujf), where I = 1,... ,no — 1. The recursion also 
shows that pnfiu^) G TLd{p){ujf). Thus G{u) G 7-frf(p)('u+)[[ wq]], and the proposition is 
proved. □ 

Remarks, i) For a solution G{u) of a g-functional equation satisfying G(0) = 0, we will 
always assume d{G) = d{P) in the following. 

ii) For the solution of a simple g-functional equation, an explicit expression may be given. 
This is e.g. the case for some g-difference equations eqn. fll.2p . with P(x, y) linear in y, see 
e.g. [H], or with P{x,y) quadratic in y, see [H]. 

Hi) It follows from G{u) G Hd{p){u+)[[uo\] that G(uo) G C[[uo]]- For derivatives of G{u), 
which are also elements of Hd{p){u+)[[uQ\], the same conclusion holds. Thus, the derivatives 




1 


G.(n) 


U=ZLo 


( 2 . 9 ) 


are formal power series, i.e., g^iuo) £ C[[mo]]- They are called factorial moment generating 
functions, for reasons to be explained in section [71 


3 Factorial moment generating functions 

Due to the following proposition, the factorial moment generating functions can be com¬ 
puted recursively, by successively differentiating the g-functional equation. 
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Proposition 3.1. Let a q-functional equation eqn. fl2.5p with solution G{u) as in Proposi- 
tion \2.5\ be given. Consider the derivative of order v ^ 0 of the q-functional equation, eval¬ 
uated at u = Uq. It is linear in gi,{uo). Its r.h.s. is a polynomial in '■ 0 ^ /j, P u}, 

with coefficients in C[[mo, 5 'o(wo)]]- 

In order to prove Proposition 13.11 we analyse partial derivatives of eqn. (12.51) . To this 
end, we employ a generalization of Faa di Brnno’s formula [15], adapted to our situation. 
The following lemma will also be used for exponent estimates in the next section. For 
k e K}, let fkiu) e C[[it]], and dehne f{u) = (/i(w),..., /a-('u)). For p G 

we use the notation f{u)^ = ((/i)p('u), ..., {fK)ti{u)) for the derivative of f{u) of order 

p. 

Lemma 3.2. (cf. [T^ Thm. 2.1]) Let a q-functional equation eqn. fl2.5p with solution G{u) 
as in Pro'position \2.5[ be given. Its derivative of order u ^ 0 satisfies 


G^{u) 


H s 

'^=1 PsG,x.) j=i 




(3.1) 


where the vectors A and hj are 1 + M + N-dimensional, the vectors puj are 1 + M- 
dimensional, for j G {1 ,..., s}, and the summation ranges over 


p,{v, X) = {(ki, ..., k,; Pi, ..., /xJ : |Ki| > 0, 

S S 

0 < Pi < ... < p^, Hi = \ and 

i=l i=l 


In the above equation, the total order < is defined as p<i/ if either |p| < \i>\, or z/|p| = \i>\, 
there exists an index fc G {0,..., M} such that pk < Vk and fVi = Vi for z = 0, ...,/c — 1. □ 

Remark. If M = 1, we have (po^Fi) <1 if and only if (pi,po) “< The 

analogous statement for higher values of M is not true. 


Proof of Proposition [Xil For given tx 7 ^ 0, we analyse the values p^ appearing in Lemma 
13.21 We show that |p^| > \i'\ for some value j , where 1 < j < s < |jz|, implies s = 1 and 
Pi = I/. Together with Lemma [2)3l this implies that pL P u. The linearity will follow from 
an explicit calculation of the term containing Pi = u. 

Assume that \pLj\ > \i'\ for some j, where 1 < j < s < \u\. The explicit form of 
A) in Lemma 13)^ states that u = This implies that \i'\ = 1^*1 ImJ- 

According to the assumption, this leads to s = 1, |k;i| = 1 and |pi| = \i'\. This, in turn, 
implies that Pi = v. 

Clearly, the r.h.s. of eqn. (13.11) . when specialised to it = tto, is a polynomial in gfj,{uo) 
for 0 7 ^ p ^ iz, with coefficients in C[[mo 5 fi'o(wo)]]- To show linearity in Gjy(tto), we note 
that the possible values of tvi are tti = ej, where j G {1,..., 1 + M + A^}. The sum of the 
terms with |pi| = \h>\ in the r.h.s. of eqn. fl3.ip is given by 

l+M+N 

i=i 
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We now extract terms containing Gi,{u) from this expression, using Lemma 1?^ and group 
them to the l.h.s. of eqn. fl3.ip . In the resulting equation, the l.h.s. L{u), when specialised 
to n = Uo, is given by 


L{uo) = (^1 - E ^(^ 0 ’ <^K)) j G^iuo). ( 3 . 2 ) 

Due to the assumptions on P{u,y), the prefactor of G,^{uo) in the above equation is 
not identically vanishing. Thus, G^{uq) is contained linearly in eqn. fl3.ip . specialised to 

U = Uq. □ 


4 Analytic generating functions 

Let Q{u, y) = P{u, 1,..., 1, y,..., y). The power series G(wo) satishes the equation 

G{uo) = Q{uo,G{uq)). (4.1) 


In the following, we specialise the class of g-functional equations. We are interested in the 
case where G(tto) is analytic at Uq = 0, with a square root as dominant singularity. This 
situation is generic for combinatorial constructions, see [HI Thm. 10.6], [161 Prop. 1], or [241 
Sec. 7.4]. Throughout the remainder of the article, we employ the following assumption. 


Assumption 4.1. Let a q-functional equation (12.5p as in Definition 12.41 be given. Let 
numbers r,s such that 0 < r, s < oo be given. Denote by Hr,d,s{u,y) the ring of power 
series in (u, y) = (mq, u+, y), which are convergent if |mo| < r, if \uk\ < d for k = 1,..., M, 
and if \yk\ < s for fc = l,...,iV. In addition to the assumptions in Definition \2.4\ we 
assume the following properties. 


i) P{u,y) G 'Hr,d,s{u,y), and its Taylor coefficients are non-negative real numbers, 
when expanded about {u, y) = (0, 0). For j = 1,..., N, each coordinate series of the 
q-shift has non-negative coefficients only, if expanded about u = 0. 


a) Let Q{u,y) = P{u,l,... ,l,y,... ,y). There exist numbers {uc,yc) satisfying 0 < 
Uc < r and 0 < yc < s, such that 


Qip'c) Vc) Vc) 


B : = 


1 ^ 
2 dy'^ 


{uc,y) 


> 0 , 

y=yc 


dQ 


= 1 , 


y=yc 


r ( 


> 0 . 


(4.2) 


Hi) The solution G{u) of the q-functional equation as in Proposition \2. ,51 has the property 
that G(uo) = J2n>iPn'^o aperiodic, i.e., there exist indices 1 < i < j < k such 
that PiPjPk 7 ^ 0, while gcd(j — i,k — i) = 1. 
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Remarks, i) When restricting to n = no, the above assumptions (together with the 
assumptions in Definition 12.41) reduce to the setup for implicitly defined power series which 
is usual in enumerative combinatorics, see [411 Thm. 10.6], [THl Prop. 1], and [24l Sec. 7.4]. 
ii) Assumption 14.11 i) implies that the coefficients pn in G{u) = J2n>oP'n^^ non¬ 
negative. For combinatorial constructions, such positivity assumptions are common. How¬ 
ever, systems of functional equations, arising from a combinatorial construction with pos¬ 
itive coefficients, might be reduced to an equation of the above form having negative coef¬ 
ficients. In that situation, other types of singularity might appear. This is, e.g., the case 
for discrete meanders [39]. The methods used in this paper can be adapted to treat such 
cases. P{u,y) G Hr,d,s{u,y) implies that, for each {u,y) G x D^x , the function 
P{u,y) has a convergent series expansion about {u,y). 

Hi) Assumption 14.11 ii) implies that the dominant singularity of G(uq) is a square root. 
This type of singularity is generic for functional equations with positive coefficients. If 
P{u,y) is a polynomial in u and y^ it can be shown that Assumption 14.11 u) follows from 
Assumption 14.11 i), if Q{u,y) is of degree > 2 in y and if Q{u,0) ^ 0. By the closure 
properties of algebraic functions, it then follows, with an adaption of Proposition 13.11 that 
all factorial moment generating functions are algebraic. 

iv) Assumption 14.II Hi) ensures that G(ito) has exactly one singularity on its circle of con¬ 
vergence. The case of periodic G{uq) may be treated by a slight extension of this setup. 

We investigate analytic properties of go{uo). A function f{u) is called A-regular [20] if 
it is analytic in the indented disc A = A{uc,ri,(j)) = {u : \u\ < Uc + p, |arg(M — Uc)| > 0} for 
some real numbers Uc > 0, p > 0 and 0, where 0 < 0 < n/2. Note that Uc ^ A, where we 
employ the convention arg(O) = 0. The set of A-regular functions is closed under addition, 
multiplication, differentiation, and integration. Moreover, if f{u) 7 ^ 0 in A, then l/f{u) 
exists in A and is A-regular. The following proposition is a straightforward extension of a 
well-known result, see e.g. [HI Thm. 10.6], [161 Prop. 1], or [2ll Sec. 7.4]. 


Proposition 4.2 (cf. [HI UHl 124]). Given Assumption \4-l\ 
analytic at Uq = 0 , with radius of convergence 0 < Uc < 00 . 

A-regular, with a sguare root singularity at Uq = Uc, and a local Puiseux expansion 


the power series go{uo) is 
Its analytic continuation is 


goiuo) = goiuc) + ^ fo,i{uc - Uq) (4.3) 

1=0 


where 70 = -1/2, /o,o = -^/G/B, and go{uc) = lim^^^^- go{uo) < 00 . □ 

Remark. The coefficients /o,z, also called amplitudes, can be computed recursively from 
the functional equation eqn. (14. Ih . by inserting the representation eqn. (14. 3 h into eqn. (14. ip 
and then expanding the resulting equation in s = This technique will be ex¬ 

ploited below, when using the method of dominant balance. 

The properties of G{uq) carry over to the factorial moment generating functions ^'^(mo)- 
We will first analyse the general form of the factorial moment generating functions, and 
later provide explicit values for exponents and leading amplitudes. 
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Proposition 4.3. Let Assumption \4-l\ he satisfied. For i/ ^ 0 arbitrary, the factorial 
moment generating function g^{uo) is analytic at Uq = 0 , with radius of convergence 0 < 
Uc < oo. Its analytic continuation is A-regular. It has a local Puiseux expansion 

OO 

gu^uo) = ^ fu,i{uc - (4.4) 

1=0 

with non-vanishing leading amplitude fi^^o ^ 0 and exponent G 

Proof. We prove the proposition by induction on w.r.t. the total order -< in Dehnition l2.2[ 
Note that the proof of Proposition 13.11 yields 

gu{uo)Qi{uo, go{uo)) = Qu{uo, go{uo), (4.5) 

where Qi{u,y) = 1 — ^{u,y) is analytic for |m| < r, \y\ < s. The function Q^{u,y,y) 
is a polynomial in y and analytic for |m| < r, \y\ < s. Due to the closure properties of 
A-regular and analytic functions, both Qi(mo, 5 'o(mo)) and Qm(mo, 5'o(mo), 
are A-regular. Due to Proposition 14.21 we have Qi(mo; 5 'o(^o)) 7 ^ 0 in A, thus its inverse is 
A-regular. This implies that guiuo) is A-regular. 

Due to Proposition 14.21 the series h{s) = yc + Yl'^o is holomorphic at s = 0 

and equals goiuo) about uq = Uc, if s = y/uc — Uq. We show that eqn. fl4.5l) leads to local 
expansions of 5 '^(mo) in terms of s = s/uc — uq, which are meromorphic at s = 0. Consider 
hrst the term Qi{uo, go{uo)). The function Qi{u,y) is holomorphic at {u,y) = {uc,yc)- 
Thus, as composition of holomorphic functions, Qi(s) := Qi{uc — s^, h{s)) is holomorphic 
at s = 0. We have Qi(0) = 0 and Q)(0) = —^^{uc,yc)h'{0) 7 ^ 0. Thus, its inverse 1/Qi{s) 
is meromorphic at s = 0, with a simple pole. Consider hnally the function Q,^{u,y,y). It 
is a polynomial in y and analytic at {u,y) = {uc,yc)- Thus, after inserting the expansions 
of g^{uo), where we use Proposition 14.21 and the induction hypothesis, we conclude that 
it is meromorphic in s at s = 0. It follows that gi,{uo) is meromorphic in s at s = 0. 
Since G{u) is non-negative and goiutf) does not vanish identically, g^iuo) does not vanish 
identically. We thus have local expansions f|4.4p of g^iuo). □ 

Remarks, i) It was argued in the preceding proof that 1 — G{uq)) = Qi{s) 

satishes Qi(s) ~ A^Uc — uq as uq ^ u~, for some constant A ^ 0, see also eqn. (13.21) . 
This will be used in the proof of the following proposition. 

a) The exponents and, in principle, all amplitudes fi,^i in the Puiseux expansion (14.41) of 
g^{uo) can be computed recursively from the functional equation eqn. fl4.1l) . compare the 
argument above for the case u = 0. Our aim in this section is to determine the exponent 
7 iy and the amplitudes fi^^. The method of Section [5] will allow to also obtain information 
about the amplitudes fi^^i for higher values of I, with reasonable effort. 

We will hrst give an estimate of 
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Proposition 4.4. Let Assumption \4-l\ be satisfied. The exponent 7 iy in the Puiseux ex¬ 
pansion of g^iuo) eqns. fl4.3D . fl4.4p . satisfies the estimate 

M . .. 

7*. < “2 + y (4-6) 

Remark. Under mild additional assumptions, the above estimate is sharp, i.e., we have 
7i^ = -^+YliLo (i + 1 ) z/j. This follows from Proposition 14.51 where we show that the am¬ 
plitudes fjyfi, for the special choice of 7 ^ as above, have non-zero values. Similar estimates 
can be obtained in situations different from a square root singularity. The proof below can 
be adapted to treat these situations. 

Proof. Set 7 jy = — r -|- where r = 1/2, and rj = (1 -|- i/2) for i = 0,..., M. We 

have rfc_|_i > for k = 0,...,M — 1. For u = 0 and u = Cq, we get from Proposition 
14.21 that 7 jy = 7 ^. We prove the proposition by induction on u w.r.t. the total order -< of 
Dehnition 12.21 using Proposition 14.31 For the remaining induction step, it suffices to show 
that 7 ^ < 7 jy. 

Assume that 7 ^ < 7 ^ holds for p ^ i/. We group all terms in eqn. fl3.ip . evaluated at 
u = Uq, which contain to the l.h.s. In the resulting equation, its l.h.s. L{s) satishes 

asymptotically, with a non-zero constant, L{s) ~ — as Uq ^ uf, see the 

remark following the proof of Proposition 14.31 By the reasoning in the proof of Proposition 
14.31 the r.h.s. R{s) of the equation satishes asymptotically, with cr a nonzero constant, 
R{s) ~ cr{uc — as uq —>• uf. Using the induction assumption, the exponent 7 can 

be estimated by 

^ / 1 ^ ^ , |A| A (|A|-1) ^ 

7 < ^ “2 + m =-^-+ 7 .- 

i=l i=l V j=0 / j=0 

In the above expression, we used Lemma [2.31 and the properties of Ps{i^, A) in Lemma [3.21 
Since L{s) equals R{s), we have 7 ^ — 1/2 = 7 . Thus, the estimate ji, < 7 ^ is satished for 
|A| > 2. Let us analyse the terms contributing to |A| = 1 in Lemma [3.21 We have s = 1, 
|ki| = 1 and = u. Now, use Lemma 12.31 If in eqn. fl 2 .ip we have |/x| = \i/\ — I for 
some I G {1,..., \v\}, it follows that 7 ,^ < 7 i/ — ^ niino<fc<M{’"fe} = lu — I- Thus, terms 
with exponent larger than 7 ^ — 1/2 must have / = 0. If |^| = \i'\ but /x 7 ^ jz, we have 
7/j, < 7iy — \'fj — 'f'k\ for some indices j,k G {0,...,M} satisfying j 7 ^ k. This means 
that the exponent estimate 7 ^ — 1/2 < — 1/2 -|- 7 ^ is satished as long as \rj — r^l > 1/2 
for j,k G {0,..., M} and j 7 ^ k. This condition is satished for the particular choice of 
exponents r* = (1 -|- i/2) of the proposition. Thus, the proposition has been proved. □ 


The reasoning in the proof of Proposition 14.41 can be rehned in order to obtain a recur¬ 
sion for the amplitudes fi^^. For vectors x = {xi,... ,xk) G and y = (?/i,..., G 
M^, we write x < y if Xi < yi for i = 1,..., K. In the proof, we will use the “large 
oh” symbol; For functions f{x) and h{x), we write f{x) = 0{h[x)) as a; —> Xq in a do¬ 
main D, if there exists a positive constant C and a neighbourhood N{xq) of Xq such that 
\f{x)\ < C\h{x)\ for all x G N{xo) D D. 
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Proposition 4.5. Let Assumption \4- 1\ he satisfied and fix u G Nq^^. Then, the amplitudes 
fv,o = fu of the Puiseux expansion eqn. fl4.4l) are, if u ^ 0 and u 7 ^ Bq, determined by the 
recursion 


M-l 


fv — Lolu-eifu-ei ^ ^ P fyf-ej+i+ej ^ ^ fpfu-pj (4-7) 


2=1 


PT^O.p^i/ 

0 <p<i/ 


with boundary conditions fo = —\JC/B < 0, /ep = —/o/2 > 0, and fiy = 0if Vj < 0 for 
some j G {1,..., M}. We have 


M 




1=0 


Vi 


For i G {0,..., M — 1}, we have /i* = —Ai/{2Bfo) > 0, and the number Ai >0 is given by 

^ 8'iP'^ BP 


i=i 


where Uc = {uc, 1,..., 1). If u Q, the amplitudes satisfy f^, > 0. This inequality is strict, 
if Ai > 0 for i = 0,..., M — 1. 

Proof. The amplitude /o has been determined in Proposition 14.21 Assume that u ^ Q and 
V 7 ^ Bq. We group all terms in eqn. fl3.ip . evaluated at tt = uq, which contain g^iuo), to 
the l.h.s. Using eqns. (13.21) and (12.91) . the l.h.s. L{uq) of the resulting equation is given by 

( dU 

Here, we used the identity 

(no) + P^fiuQ,G{uQ)), 


dG 

dtin 


N 


K) = I 5^^(no,G(no)) 


dG 

dtin 



which is obtained from differentiating eqn. (14.11) . By the reasoning in the proof of Propo¬ 
sition 14.41 the r.h.s. R{uq) of the resulting equation satishes asymptotically, with cr a 
nonzero constant, R{uq) ~ cr{uc — as uq —>• u~ . We collect all terms with 

exponents 7 iy — 1/2. Due to the proof of Proposition 14.41 they arise from terms where 
|A| = 1 or |A| = 2 in eqn. fl3.ip . An explicit analysis using Lemma [2.31 whose details we 
omit, leads to the following expressions. The contribution arising from terms where | A| = 2 
is given by 

if ^ d^P \ 

2 9p{uo)gu-p{uo). 

0 <p<i/ 
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The contribution from terms where |A| = 1 is given by 


^ f)P 

G{uo)) + ei)!-^(no)^^+e,-e,+iM- 

j=i 1=0 *+^ 

Omitting arguments and normalising the l.h.s., we arrive at the equation 

y~^ d'^P 


9u 


dyjdvk r^i (\ , G' dP OVq^ , 


2P, 


ei % dUi 


.(J) 




as Uo —> , where ()' denotes differentiation w.r.t. Uq- We have G' = ge-^ and f^i = 

-/o/2 > 0. The implied recursion for the amplitudes /j^ o is given by 

r _ ^ dyjdyk f f(i\ fS^' f f ^ fo ^ f 

“ 2Pei V 2 J \2^jpy-p) 2Pei ^ % du^ 

After rewriting prefactors, using eqn. fl4.2p and Proposition 14.31 we arrive at eqn. fl4.7|) . 

We show strict positivity of the amplitudes, if Aj > 0 for z = 0,..., M — 1. For li^l >2, 
all prefactors of in eqn. fl4.7p are non-negative. For the hrst term, this follows from the 
estimate 

/ 1 3 Az + 2 \ , , 

7„-e.= p-j-j + E—l-.j >-2 +k|>0 . 

Moreover, the last sum in eqn. 04.71) is over a nonempty set of indices p and has a strictly 
positive prefactor. For \h>\ = 1 and u ^ ei, the first term in eqn. 04.7p is zero, whereas the 
prefactors of in the second term are strictly positive. Note hnally that /e^ = Polofo > 
0, as is readily inferred from eqn. 04.71) . This leads, by induction, to strictly positive 
amplitudes The same argument shows that > 0 for ly ^ 0. □ 


5 Method of dominant balance 

We now discuss an altervative method to obtain the recursion for the amplitudes fi, in 
eqn. 04.7p . if z/q = 0. It is based on a generating function approach, and generally easier to 
apply than the combinatorial approach in the proof of Proposition 14.51 which was based on 
an application of Faa di Bruno’s formula. It also allows to analyse corrections to the leading 
singular behaviour of the factorial moment generating functions. Introduce parameters 
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6i = 1 — Ui, where i G {1,..., M}, and set <5 = ((5i,..., 6m)- From now on, we consider 
factorial moment generating functions g^{uo) for = 0 only. For k = {ki,, kM) G , 
set gk{uo) = g{o,k){uo), fk{uo) = f{o,k){uo), and 7fc(Mo) = 7(o,fc)(Mo)- 

Proposition 5.1. Assume that the q-functional equation fl2.5p has the solution G{uo, u^) G 
Hd{p){u^)[[uo\] such that (7(0,^+) = 0. Then G{uo,8) := G(mo, 1 — (5i,..., 1 — 5m) £ 
C[[mo,<5]] is a formal power series, given by 

fe >0 

where gk{uo) = (7(o,fc)(wo) ^ C[[mo]] are the factorial moment generating functions fl2.9l) . 

Proof. Proposition l2.5l states that G(u) G 7-fd(p)('it+)[[Mo]]- Thus (^(mo, <5) G 7irf(p)_i(5)[[Mo]] C 
C[[5]][[mo]] = C[[mo,< 5]]. We thus have 

G{uo,6) = Y,hk{uo)6\ 

k>0 


for some hk{uo) G C[[mo]]. By Taylor’s formula, we have 


hk{uo) — ■^G(o,fc)(wo, <5) 


5=0 


(-1)^ 

k\ 


'G(o,k){ 


u 


= (-i)VM, 


u=uo 


and the proposition is proved. 


□ 


Let C((s)) denote the held of formal Laurent series /(s) = J2i>io where Iq E T> and 
fiEC for I > Iq. Employing the Puiseux expansions of the factorial moment generating 
functions, we have an alternative representation of the generating function G{u). 


Proposition 5.2. Let Assumption \4.1\ he satisfied. Replace the coefficients gk{uo) of 
G{uo, d) G C[[mo, (5]] in Proposition I5.il by their Puiseux expansions of Propositions ff-Sf 
and \4.A o-nd denote the resulting series by G{s, 6 ) G C((s))[[5]]. It is explicitly given by 


G{s, 6 ) — G{uc) + ^ ( y^(- 


fc >0 \l=0 




where jk = where the numbers fk,i = f{o,k),i are defined in eqn. (14. 3 h 

and in eqn. fl4.4p . and where Uc = {uc, 

Then, the rescaled series G(u(s,e)) = G(s, Cis^, € 28 ^,..., eMS^~^^) G C[[s, e]] is a for¬ 
mal power series, 

G(u(s, e)) = G(uc) + sF (s, e), (5.1) 

where F{s, e) G C[[s, e]] is given by 

00 

F{s,e) = Fi{e) = ■ (5.2) 

1=0 k >0 
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Proof. Due to Propositions 14.21 and 14.31 we have gk{uo) G C((s)), where s = y/Uc — uq- 
The explicit form fl5.ll) follows immediately from the Puiseux expansions in Propositions 
14.21 and l4.3[ together with the exponent estimate in Proposition 14.41 □ 


Remarks, i) Equations like eqn. fl5.ll) appear in statistical mechanics as a so-called scaling 
Ansatz, being an assumption on the behaviour of a generating function near a multicriti- 
cal point singularity [26] , see also [121 E2] • Its validity has been proved only in a limited 
number of examples. Here, we employ the different framework of formal power series. 
a) Proposition 15.21 suggests an alternative strategy to compute the singular behaviour 
of the factorial moment generating functions. We consider the functional equation for 
F{s,e) G C[[s, e]], which is induced by the g-functional equation for G{u). Writing 
F(s, e) = where Pz(e) G C[[e]], then leads to a partial differential equa¬ 

tion for Fi{e), upon expanding the induced functional equation in powers of s. Note that 
this algorithm, which is computationally involved, can be easily implemented in a com¬ 
puter algebra system. The above method is called in statistical mechanics the method of 
dominant balance, see [441IT6] . It can be applied in our framework, if an exponent bound 
like that of Proposition 14.41 is known. Such bounds are generally easier to obtain than 
explicit recursions for amplitudes, see the above proofs. 


Proposition 5.3. Let Assumption \4.1 
{uo{s, e),..., mm(s, e)), where 


be satisfied. 


Define the power series u{s, e) = 


uo{s,e) = Uc-s‘^, Mi(s, e) = 1 - (z = l,...,M). 


For a q-shift v, define the induced q-shift of s and the induced q-shift e^, = ..., eM,v) 

ofe by 


S'u 


= S^{s, e) ■■ 


\/mc - Voiu{s, e)), 
1 - Vi{u{s, e)) 
s^,(s, e)*+2 


(z = l,...,M). 


We then have Sjj G C[[s, e]] and G C[[s, e]] for i = 1,..., M. The functional equation 
for F{s,e) G C[[s, e]] of Proposition IFTB. induced by eqn. fl2.5p . is given by 


G{u{s, e)) = P{u{s, e), H{u{s, e))). (5.3) 

In the above equation, G{u{s, e)) is given by eqn. fl5.1l) . and the power series e)) 

are given by 

e)) = G{ufi) + e^o)) (1 < j < N). 


Proof. This follows from direct computation. Note hrst that the Taylor expansions of the 
power series Vk{u) about u = Uq are of the form 

Vk{u) =Vk{uo)+ ^^-^^Vk,i{uo){l-Uiy^ ■ ...■ (fc = 0,...,M). (5.4) 

1^0 

l>0 
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Now insert the parametrisations Uk{s, e) G C[s, e], where k = 0,, M. We get Vk{u{s, e)) G 
C[[s, e]]. Due to the g-shift properties, we have 

Vo{u{s, e)) =Uc- - ^{uc)eis^ + 0 (/), 

uU\ 

Vk{u{s, e)) = 1 - {u^)ek+is'^""^ + C>(s^+^) (A; = 1,..., M - 1), 

OUk+i 

vm{,u{s, e)) = 1 - + C>(s^+^). 


We thus have Uo('*^('S, e)) = Uc — — s^Ro{s,e) and Vk{u{s,e)) = 1 — s^~^^Rk{s, e) for 

k = 1,. .., M, where Ri{s, e) G C[[s, e]] for i = 0,..., M. For s^, we get 

Sv = \/uc- Vo{u{s, e)) = sa/1 - sRo{s, e). 

We thus have G C[[s, e]] and s„(0, e) = 0. We get 

e)) Ri{s,e) 




^2 + 2 


— sR{s, e) 


i+2 




It follows that e^,^, G C[[s, e]], where i = 1,..., M. Furthermore, we get from eqn. (15.dh 

ei,«(0, e) = i?i(0, e) = e^. 


such that € 1 ,( 0 , e) = e. This ensures that eqn. fl5.3p is well-dehned for F{s, e) G C[[s, e]]. □ 

Proposition 5.4. Let Assumption \4.1\ he satisfied. The formal power series Fo{^) G C[[e]] 
satisfies the singular, quasi-linear partial differential equation of first order 


N 


dP 


1 dv, 


U) 


M 


d 


(n„G(n.)) + ) +BF^{e)-C = 0, (5.5) 


2=1 


where the numbers B, C are defined in eqn. (14.2h . and the formal power series hl^\e) are 
given by 


tdHe) 

Rho 


dyiL 

dui+i 


Mei+i 


M + 2 dv\f^ 
2 dui 


i + 2 dv^'^ 
2 dui 


{ufieiCi 


{uc)eieM- 




(5.6) 


Proof. Using the expansions in the proof of Proposition 15.31 it is readily inferred that the 
power series s„ and e.„ are, to leading orders in s, given by 


1 

SijO) = s + - „ ° {uc)eis‘^ + (P(s^), 

2 oui 

^i,vU) = + F hl^\e)s + 0(3“^) (z = 1,..., M — 1), 
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with h^P{e) as defined in eqn. (I5.6p . Using this resnlt, we compnte the expansion of 
e)) np to order This yields 

e)) = G{uc) + s^,o■)-^(s^,o■), e^O)) 

= G{uc) + FQ{e)s + f -Fi(e) + - ('Uc)eiFo(e) + hp'*(e)—Fo(e) J + O(s^). 

Now expand the fnnctional eqnation fl5.3p to leading orders in s. Terms of order vanish 
dne to eqn. fl4.ll) . evalnated at Uq = Uc- Terms of order vanish dne the condition 
§p(^c, G{Uc)) = 1. Terms of order lead to the partial differential equation given 
above. □ 


Remarks, i) For M > 2, it is an open question whether closed form solutions for Fo{e) 
exist. See Ha for a discussion of the cases M = 1 and M = 2. 

ii) The above method can also be used to derive partial differential equations characterising 
the generating functions Fi{s, e) of the amplitudes fk,i for I > 0. These equations arise in 
the expansion of the g-functional equation in s at order 1 + 2, see [16] for examples where 
M = 1. 

The above theorem leads to an alternative derivation of the recursion eqn. (14.71) of 
Proposition 14.51 in the case z/q = 0. 


Alternative proof of eqn. fl4.7p . We set Fo(e) = K{—e) + /o,o and rewrite eqn. fl5.5l) in the 
form 


rrt ^ i + 2. dK{e) 

- L —'‘"'“‘“ay 

2 = 1 


^ dK{e) 
2=1 


dei 


^ei{K{e) + /o,o) - 


-Kief, (5.7) 


where K{0) = 0, and the constants fii are, for i G {0,...,M — 1}, given by /ij = 
—Ai/{2Bfo^o), with Ai defined in eqn. (14.8p . This leads to the recursion eqn. (14.71) for 
the coeffients /jy in Proposition 14.51 if we set u = (0, k). □ 


6 Growth of amplitudes 


We are interested in the growth of the coefficents fi^^, which appear in Proposition 14.51 
and in Proposition 15.41 in the case vq = 0. To this end, we study properties of the partial 
differential equation of the associated generating function Fo(e) of eqn. fl5.2|) . given by 


f’o(c) = 


fc >0 


Proposition 6.1. For a q-functional equation, let Assumption \4.1\ he satisfied. There exist 
positive real numbers D, Ri ,..., Rm, such that 

|/fc,ol < D{k^ + • • • + kM)\{Rif^ ■... ■ {Ruf^ 


for all k > 0. 
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For the proof of the proposition, we apply the technique of majorising series. A formal 
power series g = majorises a formal power series h = if \gk\ < \k\\hk 

for all fc > 0. We then write g <^h. We have the following relations. 

Lemma 6.2. Let g ^h. Then 

212 1 dg ^ 

g h , xig -C Xih, XiXi -— -C xih, 

UXi 

for i = 1,..., M and j = 1,..., M — 1. 

Sketch of proof. These relations are checked by direct computation. If i = 2,..., M, we 
have XiXi^, = Thus \kigk-ei\ < (1^1 - ^Y-kikk-e^ < |fc|!hfc-ei, and the 

statement follows. The remaining assertations are proved similarly. □ 

Proof of Proposition [Ail Starting with K{e), dehned in eqn. fl5.7p . we introduce the ma- 
jorant equation 

T(e) = po 2 ^ ‘^iT(e) + ^ei(L(e) + |/o,o|) + ^j-^-^T(e)^+ ^ p-iCi+i ^ ■ (6.2) 

The last equation uniquely dehnes a power series with non-negative coefficients and pos¬ 
itive radius of convergence, satisfying L{0) = 0. It belongs to a class of singular partial 
differential equations with regular solutions, which is discussed in [2S1 Thm. 2.8.2.1 and 
Sec. 2.9.5]. To prove that K{e) T(e), denote the r.h.s. of eqn. fl5.7p by RK and the 
r.h.s. of eqn. fl6.2p by RL. By construction, g ^ h implies Rg Rh. We use an iteration 
argument. Set Kq = Lq = 0. Clearly Kq -C Lq. Dehne Kn = RKn-i and = RLn-i 
for n G N. We have Kn -C for n G Mq, due to Lemma [621 Thus K{e) -C L{e) for the 
formal solutions K{e) of eqn. fl5.7p and L(e) of eqn. 06.21) . Since L(e) is regular, we have 
the estimate 

|/fc,o| < \k\\[e’^]L{e) < DlkllR’^ 

for some real constants D > 0 and Ri > 0, where i = 1,, M. □ 

7 Moments and limit distributions 

In the following, we give a probabilistic interpretation of the obtained results. This gener¬ 
alises the discussion of Dyck paths in the introduction, before Proposition 11.41 to the case 
of a general g-functional equation, and will prove part i) of Theorem 11.51 For technical 
reasons (Lemma 17.4p . we will not use random variables below, but rather argue with the 
associated probability measures. 

For a g-functional equation, let Assumption 14.11 be satished. Then, the coefficients pn 
in a solution G{u) = 'YhnyoVn'uP, such that G(0) = 0, are non-negative. Assume that the 
numbers Ai of Proposition 14.51 satisfy Aj > 0 for i = 0,..., M — 1. We then have fkp > 0 


dg dh , . 
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for all fc > 0 and fc 7 ^ 0. This implies that nMP^o,ni,...,nM > 0 almost all uq, 

see eqn. fl7.2p below. Fix iVo G N such that strict positivity holds for all Hq > Nq. For 
no > No, we dehne discrete Borel probability measures Jlno by 

~ _ \ ^ Pno,ni,...,nM r 

Pno — 2 ^ ^ 0(ni,...,nM)) 


compare eqn. fll.lSp in the introduction. Their corresponding moments are, for k G 
given by 


mfc(no) 


E 




n 


ki 


Pno,ni,...,nM 


E 




Pno 




(7.1) 


We are interested in the asymptotic behaviour of the moments fhk{no), as no tends to 
inhnity. This will be obtained by an analysis of the coefficients of the factorial moment 
generating functions, which relies on a transfer lemma [22l Thm. 1]. 


Lemma 7.1 (Transfer lemma [22]). Suppose that F{z) = En>o ® singularity at 

z = Zc and is A-regular, i.e., it is analytic in the domain 


A = {z ■.\z\ < Zc + f], \aTg{z - Zc)\ > 0} 
for some rj > 0 and 0 < (f < 7r/2. Assume that, as z ^ Zc in A, 

F{z) = 0{{z,-zr) 

for some real a. Then, the n-th Taylor coefficient fn of F{z) satisfies 


u = Oizf^n-^-^) {n ^ 00 ). 


□ 

The following lemma characterises the asymptotic behaviour of the moments mfc(no), 
as no tends to inhnity. 

Lemma 7.2. Let Assumption lf. 1\ be satisfied. Assume that the numbers Ai of Proposition 
\4.5\ satisfy Aj > 0 for i = 0,..., M — 1. Then the moments rhk{no) egn. fl7.ll) are well- 
defined for almost all no- They are for k G asymptotically given by 


mfc(no) = 


k\ 


r(7o) 


r(7fc) 


+ Oinl 


(no 


00 


where F(z) denotes the Gamma function, and where the numbers f^p and defined 

in Proposition\4.5\ 
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Proof. The functions gk{uo) are A-regular due to Proposition I4.3[ We infer from eqn. fl4.4p 
that gk{uo) = fk,oiuc — + 0 {{uc — as Uq —>• u~, where 70 = — 1/2 and 

7 fc > 0 otherwise. The coefficient asymptotics of {uc — and the transfer lemma yield 

[uo°]9k{uo) = (1 + C>(no^^^)j (no^oo), (7.2) 

where [x^]f{x) denotes the coefficient of a:” in the power series f{x). The error term implies 
that asymptotically factorial moments coincide with ordinary moments. The numbers 
T^kipo) are thus well-dehned for no large enough, and asymptotically given by 


mfc(no) = 


'n,n+>Q'^XPno,n+ _ [uQ°]gk{uQ) 


'Yl,n+>QPno,n+ 
k\ r(7o) 
/o,ow?“^“ r( 7 fc) 


k\ 


«o J^o(mo) V J 

(^1 + C>(no j (no 


CX)). 


This concludes the proof of the lemma. 


□ 


The moments mk{no) diverge as no —> 00 , as may be inferred from Lemma [7.21 Intro¬ 
duce normalised Borel probability measures fino by 


Pno 




Pno, 






E 




Pno 






(7.3) 


where n^ = n^np ^ _ x,..., M. For k G denote the corresponding moments 

by nrfc(no). We now show that the limits 


rrik = lim mk{no) = lim 

TiQ—>00 no—^0( 


mfc(no) 


k\ 


r(7o) 


n. 


7fc-70 

0 


r(7fc) 


/fc,( 


(7.4) 


exist and dehne a unique Borel probability measure /r, with hnite moments ruk > 0 at all 
orders. This will be achieved using Levy’s continuity theorem. We hrst prove the following 
lemma. 


Lemma 7.3. Let Assumption 4-1 be satisfied. Assume that the numbers Ai of Proposition 
4 . 5\ satisfy Ai > 0 fori = 0,..., M — 1. Consider for k G the numbers mk > 0, defined 
in eqn. fl7.4l) and eqn. fl7.ip . Fort G we have 


lim 

|fc|—^OO 


nikt^ 

k\ 


0 . 


Equivalently, the exponential generating function of the numbers mk is entire. 
Proof. The limit mk eqn. (17.41) exists for k G , due to Lemma [7.21 Note that 


7fc 2 


M 

E 

2=1 


i\ 13 

IH— 1 ki > -1— Ifcl > I 

2 2 2' ^ 


\k\ 


-J- 
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Thus |fc| —>■ cx) implies 7 fc —> oo. Furthermore, since einjeY' ^ n\ < en{n/e)'^ for n G M, 
we have the estimate 


n\ 


< 


ninjey 




< 


(n + uo)! ((u + no)/e)”+’^o n^o-^ 

for u, no G N. Now hx t G We estimate 


rukt^ 


k\ 


< 


r(7o) I/mI 

r(7j /o,o«?‘-''” 

/o,o 


t^\ < 


r(7o) 


D\k\\W 


(ifci + m)! /o.o 


It'"! max(Mc, W ) 


-l'\7fe-70 


\k\^' 


2 J 


-1 


M+2 


|fc| 


The r.h.s. vanishes as \k\ —> oo, which implies the assertion. The equivalent statement is 
obvious. □ 


Our proof of claim i) in Theorem 11.51 relies on an application of Levy’s continuity 
theorem [5l Thm. 23.8], which we cite in the following lemma. 

Lemma 7.4 (Levy’s continuity theorem [5]). Forn G N, let probability measures fin on the 
Borel a-algebra ofM.^ be given. If the seguence {/in}neN of their characteristic functions 
fin converges pointwise to a complex function (f which is continuous at the origin, then cf is 
the characteristic function of a uniguely determined Borel probability measure fi. Moreover, 
{h'nfneN converges to fi weakly. □ 

Proof of claim i) in Theorem \1.5[ Lemma [7.21 implies that the Borel probability measures 
fing eqn. fl7.3p are well-dehned for almost all uq G M. The estimate in Lemma 17.31 implies 
uniform convergence of the sequence of Fourier transforms fino '■ —>■ C of fin^, in every 

ball of hnite radius centred at the origin. In particular, we have pointwise convergence of 
the sequence {/inojnoeN- For M = 1, the corresponding argument is given in the proof of 
m Thm 6.4.5]. It can be directly extended to arbitrary M. Since the functions fino are 
continuous, we conclude that the limit function 0 : —>■ C is continuous at the origin. 

Now apply Levy’s continuity theorem. The limit probability measure fi has moments 
eqn. fl7.4p . The claimed statement of the theorem follows, when phrasing the result in 
terms of the associated random variables. □ 


The connection to Brownian motion, which is claimed in part ii) of Theorem 11.51 will 
be established in the following section. 


8 Dyck paths revisited 

We apply our results to the example of Dyck paths of Section [L^ The power series solution 
E{uo) of eqn. fl2.6p . specialised to u = uq, has radius of convergence Uc = 1/4, with a 
square root singularity at u = Uc, and E{uf) = 1. The g-functional equation eqn. fl2.6p 
satishes Assumption 14.11 Since the random variables (Wi ..., XM,no) of 111-51) have 
the distribution of fin^, as dehned in eqn. fl7.3l) . Proposition 11.41 follows from the results of 
the previous section. 


30 



































Proof of Proposition\1.4\ The generating function E{u) of Dyck paths satishes the q- 


functional equation eqn. fl2.6l) . Assumption 14.II holds with Uc = 1/4 and Hc = 1. We have 
/o,o = —4, 7o = —1/2, fii = {i + l)/4 for z = 1,..., M — 1 and qiQ = 1/8. Thus, Propo¬ 
sition 02] yields eqn. (11.111) . The distribution of the random variables ... ,XM,no) 

of eqn. (11.51) is that of the probability measures qino in eqn. (17.dh . By part i) of Theo¬ 
rem 11.51 there exists a unique limit probability measure /z. We thus get eqn. fll.161) from 
eqn. fl7.4p . □ 


The connection between Dyck paths and Brownian excursions in Proposition 11.11 leads, 
for a general g-functional equation, to an explicit characterisation of the limit probability 
measure /z, resp. of the associated limit random variable (Ti,..., Ym)- 

Proof of claim ii) in Theorem \1.5[ For the given g-functional equation, let Fo(e) be the 
generating function of the leading amplitudes f^^. For k = 1,... ,M and dk € M, dehne 
Go{e) = pQ^eidi ,..., CMdM)- An explicit calculation using eqn. fl5.7p shows that the power 
series G'o(e) satishes the same type of differential equation as Fo(e) does, with /zq replaced 
by fTodi and /ij replaced by /Zjdi+i/dj, where z = 1,..., M — 1. Now choose the values dk, 
such that the equation for Go{e) is that of Dyck paths. Noting the relation between Dyck 
paths and Brownian excursions eqn. fll.7p . we arrive at the values Ck = /dk, for 

numbers Ck as in the claim of Theorem 11.51 □ 


9 Concluding remarks 

We hnally stress three central aspects of our approach. Firstly, the approach yields a 
universal limit distribution - loosely spoken, it appears for all models, whose underlying 
functional equation has the same singularity structure. Such a result may be compared to 
a central limit theorem in probability theory. For example, models other than Dyck paths, 
which display a square root as dominant singularity of their size generating function, are 
certain models of trees or polygons. For simply generated trees, ^-functional equations 
appear when counting by number of vertices and by internal path length [50]. Using 
the above setup, moment recursions for the parameter “sum of fc-th powers of the vertex 
distances to the root” are obtained. For polygon models, g-functional equations appear 
when counting by perimeter and area [T^ 14^ . which is, for column-convex polygons, the 
sum of the column heights. The above setup gives moment recursions for the parameter 
“sums of fc-th powers of the column heights”, in the limit of inhnite (horizontal) perimeter. 

Secondly, our approach is algorithmic - the moment recursion can hnally be deduced 
from a straightforward calculation, by the method of dominant balance. Our approach 
also allows to obtain corrections to the asymptotic behaviour, which cannot (easily) be 
deduced by other methods. 

Thirdly, the approach is flexible - it may be applied to other classes to obtain a uni¬ 
versal limit distribution, which only depends on the singularity structure of the functional 
equation. Our generating function approach is particularly suited for counting parame¬ 
ters, which decompose linearly under the cartesian product construction. Examples of 
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such models with a rational generating function appear in [16]. Examples with an inverse 
square root appear in [39]. In particular, the discrete counterparts of Brownian motion, 
Brownian bridges, and Brownian meanders can be studied, see [39l [48] . Also, universality 
questions for parameters related to left and right path lengths in trees [29], [33l Il2l HD] can 
be studied by our methods, compare the discussion in m- 

Within the framework of simply generated trees, an alternative derivation of the above 
moment recursion could be obtained with the techniques of [28] , where the different problem 
of the (generalised) Wiener index of trees was analysed. It would be interesting to consider 
how our methods can be adapted to this problem. 


Acknowledgements 

The author thanks Philippe Duchon, Philippe Flajolet and Michel Nguyen The for helpful 
discussions, and Svante Janson for comments on the manuscript. He thanks the department 
LaBRI (Bordeaux) for hospitality in autumn 2003, where parts of the problem have been 
analysed. The author thanks the referees for suggestions, which improved the presentation 
of the article. Financial support by the German Research Council (DFG) is gratefully 
acknowledged. 


References 

[1] D.J. Aldous, The continuum random tree I, Ann. Prob. 19 (1991), 1-28. 

[2] D.J. Aldous, The continuum random tree II: an overview. Stochastic Analysis, eds. 
M.T. Barlow and N.H. Bingham, Cambridge University Press (1991), 23-70. 

[3] D.J. Aldous, The continuum random tree III, Ann. Prob. 21 (1993), 248-289. 

[4] C. Banderier and P. Flajolet, Basic analytic combinatorics of directed lattice paths, 
Theoret. Comput. Sci. 281 (2002), 37-80. 

[5] H. Bauer, Probability Theory, de Gruyter Studies in Mathematics, vol. 23, de Gruyter, 
Berlin (1996). 

[6] B. Belkin, An invariance principle for conditioned recurrent random walk attracted to 
a stable law, Z. Wahrscheinlichkeitstheorie und Verw. Gebiete 21 (1972), 45-64. 

[7] P. Billingsley, Probability and Measure, 3rd ed., Wiley, New York (1995). 

[8] P. Billingsley, Convergence of Probability Measures, 2nd ed., Wiley, New York (1999). 

[9] M. Bousquet-Melon, A method for the enumeration of various classes of column-convex 
polygons, Discr. Math. 154 (1996), 1-25. 


32 


[10] M. Bousquet-Melou, Limit laws for embedded trees: applications to the integrated 
superBrownian excursion, Random Structures Algorithms 29 (2006), 475-523. 

[11] M. Bousquet-Melou and S. Janson, The density of the ISE and local limit laws for 
embedded trees, Ann. Appl. Probab. 16 (2006), 1597-1632. 

[12] J. Bouttier, P. Di Francesco and E. Quitter, Geodesic distance in planar graphs. 
Nuclear Rhys. B 663 (2003), 535-567. 

[13] N.G. de Bruijn, D.E. Knuth, and S.O. Rice, The average height of planted plane trees, 
in: Graph Theory and Computing, pp. 15-22, Academic Press, New York (1972). 

[14] K.L. Ghung, A Gourse in Probability Theory, Academic Press, New York (1974). 

[15] G.M. Gonstantine and T.H. Savits, A multivariate Faa di Bruno formula with appli¬ 
cations, Trans. Amer. Math. Soc. 348 (1996), 503-520. 

[16] M. Drmota, Systems of functional equations. Random Struct. Algorithms 10 (1997), 
103-124. 

[17] M. Drmota, Stochastic analysis of tree-like data structures, Proc. Royal Soc. Lond. A 
460 (2004), 271-307. 

[18] P. Duchon, Q-grammaires: un outil pour I’enumeration, PhD thesis, Bordeaux Uni¬ 
versity (1998). 

[19] P. Duchon, Q-grammars and wall polyominoes, Ann. Comb. 3 (1999), 311-321. 

[20] J.A. Fill, P. Flajolet and N. Kapur, Singularity analysis, Hadamard products, and 
tree recurrences, J. Comput. Appl. Math. 174 (2005), 271-313. 

[21] P. Flajolet and G. Louchard, Analytic variations on the Airy distribution, Algorithmica 
31 (2001), 361-377. 

[22] P. Flajolet and A.M. Odlyzko, Singularity analysis of generating functions, SIAM 
J. Discr. Math. 3 (1990), 216-240. 

[23] P. Flajolet, P. Poblete and A. Viola, On the analysis of linear probing hashing, Algo¬ 
rithmica 22 (1998), 37-71. 

[24] P. Flajolet and R. Sedgewick, Analytic Combinatorics, book in preparation (2007). 

[25] R. Gerard and H. Tahara, Singular Nonlinear Partial Differential Eguations, Aspects 
of Mathematics, vol. E 28, Vieweg, Braunschweig (1996). 

[26] E. J. Janse van Rensburg, The Statistical Mechanics of Interacting Walks, Polygons, 
Animals and Vesicles, Oxford University Press, Oxford (2000). 


33 



[27] S. Janson, T. Luczak and A. Rucinski, Random graphs, Wiley, New York (2000). 

[28] S. Janson, The Wiener index of simply generated random trees. Random Struct. Al¬ 
gorithms 22 (2003), 337-358. 

[29] S. Janson, Left and right pathlengths in random binary trees, Algorithmica 46 (2006), 
419-429. 

[30] S. Janson, Brownian excursion area, Wright’s constants in graph enumeration, and 
other Brownian areas, Probab. Surv. 4 (2007), 80-145. 

[31] M.J. Kearney, S.N. Majumdar and R.J. Martin, The first-paasage area for drifted 
Brownian motion and the moments of the Airy distribution, preprint (2007); 
arXiv:0706.2038. 

[32] C. Knessl and W. Szpankowski, Enumeration of binary trees and universal types. 
Discrete Math. Theor. Comput. Sci. 7 (2005), 313-400. 

[33] C. Knessl and W. Szpankowski, On the joint path length distribution in random binary 
trees. Stud. Appl. Math. 117 (2006), 109-147. 

[34] G. Louchard, Kac’s formula. Levy’s local time and Brownian excursion, J. Appl. Prob. 
21 (1984), 479-499. 

[35] G. Louchard, The Brownian excursion area: a numerical analysis, Com¬ 
put. Math. Appl. 10 (1985), 413-417; Erratum, Comput. Math. Appl. 12 (1986), 375. 

[36] J.-F. Marckert and A. Mokkadem, The depth hrst processes of Galton-Watson trees 
converge to the same Brownian excursion, Ann. Probab. 31 (2003), 1655-1678. 

[37] A. Meir and J.W. Moon, On the altitude of nodes in random trees, Canad. J. Math. 
30 (1978), 997-1015. 

[38] M. Nguyen The, Distributions de valuations sur les arbres, PhD thesis, L’Ecole Poly¬ 
technique, Paris (2003). 

[39] M. Nguyen The, Area of Brownian Motion with Generatingfunctionology, in; Discrete 
Random Walks, DRW’03, eds. G. Banderier and G. Krattenthaler, Discr. Math, and 
Theoret. Comput. Sci. Proceedings AC (2003), 229-242. 

[40] M. Nguyen The, Area and inertial moment of Dyck paths, Combin. Probab. Comput. 
13 (2004), 697-716. 

[41] A.M. Odlyzko, Asymptotic enumeration methods, in; Handbook of Combinatorics, 
vol. 2, eds. R.L. Graham, M. Grotschel and L. Lovasz, Elsevier, Amsterdam (1995), 
1063-1229. 


34 



[42] A. Panholzer, Left and right length of paths in binary trees or on a question of Knuth, 
Discr. Math, and Theoret. Comput. Sci. Proceedings AG (2006), 415-418. 

[43] T. Prellberg, Uniform g-series asymptotics for staircase polygons, J. Phys. A: 
Math. Gen. 28 (1995), 1289-1304. 

[44] T. Prellberg and R. Brak, Critical exponents from non-linear functional equations for 
partially directed cluster models, J. Stat. Phys. 78 (1995), 701-730. 

[45] T. Prellberg and A.L. Owczarek, Stacking models of vesicles and compact clusters, 
J. Stat. Phys. 80 (1995), 755-779. 

[46] C. Richard, Scaling behaviour of two-dimensional polygon models, J. Stat. Phys. 108 
(2002), 459-493. 

[47] C. Richard, Limit distributions for models of exactly solvable walks, in: Oberwol- 
fach reports 1, Report No. 22/2004, Mathematisches Forschungsinstitut Oberwolfach 
(2004), 1189-1191. 

[48] C. Richard, Staircase polygons: moments of diagonal lengths and column heights, 
J. Phys.: Gonf. Ser. 42 (2006), 239-257. 

[49] R.P. Stanley, Enumerative Gombinatorics, vol. 2, Cambridge University Press, Cam¬ 
bridge (1999). 

[50] L. Takacs, A Bernoulli excursion and its various applications, Adv. Appl. Prob. 23 
(1991), 557-585. 

[51] L. Di Vizio, J.-P. Ramis, J. Sauloy and C. Zhang, Equations aux g-differences, 
Gaz. Math. 96 (2003), 20—49. 


35 



